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Abstract. The presence of 6 Li in the atmospheres of metal-poor halo stars is usually in- 
ferred from the detection of a subtle extra depression in the red wing of the 7 Li doublet 
line at 670.8 nm. However, as pointed out recently by Cayrel et al. (2007), the intrinsic line 
asymmetry caused by convective flows in the photospheres of cool stars is almost indis- 
tinguishable from the asymmetry produced by a weak s Li blend on a (presumed) symmet- 
ric 7 Li profile. Previous determinations of the 6 Li/ 7 Li isotopic ratio based on ID model 
atmospheres, ignoring the convection-induced line asymmetry, must therefore be consid- 
ered as upper limits. By comparing synthetic ID LTE and 3D non-LTE line profiles of the 
Li 1 670.8 nm feature, we quantify the differential effect of the convective line asymmetry 
on the derived 6 Li abundance as a function of effective temperature, gravity, and metallicity. 
As expected, we find that the asymmetry effect systematically reduces the resulting 6 Li/ 7 Li 
ratios. Depending on the stellar parameters, the 3D-1D offset in 6 Li/ 7 Li ranges between 
-0.005 and -0.020. When this purely theoretical correction is taken into account for the 
Asplund et al. (2006) sample of stars, the number of significant 6 Li detections decreases 
from 9 to 5 (2<r criterion), or from 5 to 2 (3<x criterion). 

We also present preliminary results of a re-analysis of high-resolution, high S/N spectra 
of individual metal-poor turn-off stars, to see whether the second Lithium problem actually 
disappears when accounting properly for convection and non-LTE line formation in 3D 
stellar atmospheres. Out of 8 stars, HD 84937 seems to be the only significant (2<r) detection 
of 6 Li. In view of our results, the existence of a 6 Li plateau appears questionable. 

Key words, stars: abundances - stars: atmospheres - convection - line: profiles - stars: 
population II - stars: individual: (HD 74000, HD 84937, HD 140283, HD 160617, 
G271-162, G020-024, G64-12, G275-4) 



1. Introduction 

A systematic analysis of high-dispersion, high 
signal-to-noise (S/N) UVES and Keck spectra 
of about 30 bright metal-poor stars by Asplund 



et al. (2006) and Asplund & Melendez (2008) 
(henceforth A06 and A08, respectively) re- 
sulted in the detection of 6 Li (at the 2 cr level) 
in about one third of these objects. The average 
6 Li/ 7 Li isotopic ratio in the stars in which 6 Li 
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has been detected is about 4% and is very simi- 
lar in each of these stars, defining a 6 Li plateau 
at approximately log n( 6 Li)/n(H) + 12 = 0.85. 
A convincing theoretical explanation of this 
new 6 Li plateau turned out to be problematic: 
the high abundance of 6 Li at the lowest metal- 
licities cannot be explained by current models 
of galactic cosmic -ray production, even if the 
depletion of 6 Li during the pre-main-sequence 
phase is ignored (see reviews by e.g. Cayrel et 
al. 2008; Christlieb 2008; Prantzos 2010, 2012, 
and references therein). 

A possible solution of this so-called second 
Lithium problem was proposed by Cayrel et 
al. (2007): noting that the spectroscopic signa- 
ture of the presence of 6 Li in the atmospheres 
of metal-poor halo stars is just a subtle ex- 
tra depression in the red wing of the 7 Li dou- 
blet, which is difficult to measure even in spec- 
tra of the highest quality, they point out that 
the intrinsic line asymmetry caused by convec- 
tion in the photospheres of cool stars closely 
mimics the presence of 6 Li at the level of a 
few percent if interpreted in terms of ID, in- 
trinsically symmetric blend components. As a 
consequence, the 6 Li abundance derived so far 
by using ID model atmospheres, ignoring the 
convective line asymmetry, are expected to be 
systematically too high. 

We quantify the theoretical effect of the 
convection-induced line asymmetry on the re- 
sulting 6 Li abundance by fitting a given syn- 
thetic profile both with a grid of ID LTE and 
a grid of 3D non-LTE line profiles. The 1D-3D 
difference of the required 6 Li/ 7 Li isotopic ra- 
tio measures the expected systematic error in 
6 Li/ 7 Li inherent to the standard ID analysis. 
The synthetic spectra used for this differential 
comparison are based on a set of 3D hydrody- 
namical model atmospheres computed with the 
CO s BOLD code (Freytag et al. 2002, 2012) 
and a corresponding set of fully compatible 
ID mixing-length models, respectively, as out- 
lined in Sect. 2. The results of this investigation 
may be represented by an analytical approx- 
imation, giving the 3D correction of 6 Li/ 7 Li 
as a function of effective temperature, grav- 
ity, and metallicity, for a parameter range that 
covers the stars of the A06 sample (Sect. 3). 
We also discuss the potential advantages and 



disadvantages of using so called 'calibration 
lines' for fixing the residual line broadening, 
and show that the usage of such 'calibration 
lines' is potentially dangerous, because the in- 
ferred broadening parameter shows consider- 
able line-to-line variations. Depending on the 
choice of these lines, the resulting 6 Li abun- 
dance may be systematically biased (Sect. 4). 

A careful reanalysis of individual objects is 
under way to see whether the second Lithium 
problem can be resolved when accounting 
properly for convection and non-LTE line for- 
mation in 3D stellar atmospheres. In Sect. 5 we 
present a preliminary analysis of eight individ- 
ual metal-poor turn-off stars for which suffi- 
ciently high-resolution, high S/N spectra are at 
our disposal, fitting the observed Li 1 670.8 nm 
feature with both ID LTE and 3D non-LTE 
synthetic line profiles. As expected, the 3D 
analysis gives a systematically lower 6 Li/ 7 Li 
ratio by roughly 0.01 with respect to the ID 
result. In most cases, we find that the detection 
of 6 Li is not significant at the 2cr level. 

2. 3D model atmospheres and 
spectrum synthesis 

The hydrodynamical atmospheres used in the 
present study are part of the CIFIST 3D 
model atmosphere grid (Ludwig et al. 2009). 
They have been obtained from realistic nu- 
merical simulations with the C0 5 BOLD code 1 
(Freytag et al. 2002, 2012) which solves 
the time-dependent equations of compress- 
ible hydrodynamics in a constant gravity 
field together with the equations of non-local, 
frequency-dependent radiative transfer in a 
Cartesian box representative of a volume lo- 
cated at the stellar surface. The computa- 
tional domain is periodic in x and y direc- 
tion, has open top and bottom boundaries, and 
is resolved by typically 140x140x150 grid 
cells. The vertical optical depth of the box 
varies typically from log T Ross a;-7.0 (top) to 
logT Ross «+7.5 (bottom). Further information 
about the grid of models used in the present 
study is compiled in Table 1 . As indicated in 

1 http : //www . astro . uu . se/~ bf /co5bold_ 
main . html 



284 



M.Steffen et al.: 6 Li detection in metal-poor stars 



columns (2) - (4), the grid covers the range 
5900 K < r eff < 6500 K, 3.5 < \ogg < 4.5, 
-3.0 < [Fe/H] < 0.0. Radiative transfer is 
solved in 5 or 6 opacity bins (col. 5). Each of 
the models is represented by typically 20 snap- 
shots chosen from the full time sequence of the 
corresponding simulation (col. 6). 

These representative snapshots are pro- 
cessed by the non-LTE code NLTE3D that 
solves the statistical equilibrium equations for 
a 17 level lithium atom with 34 line transi- 
tions, fully taking into account the 3D thermal 
structure of the respective model atmosphere 
(but neglecting the differential Doppler shifts 
implied by the hydrodynamical velocity field). 
The photo-ionizing radiation field is computed 
at 704 frequency points between A 925 and 
32 407 A, using the opacity distribution func- 
tions of Castelli & Kurucz (2004) to allow for 
metallicity-dependent line-blanketing, includ- 
ing the HI-H + and HI-HI quasi-molecular 
absorption near A 1400 and 1600 A, respec- 
tively. Collisional ionization by neutral hydro- 
gen via the charge transfer reaction H(ls) + 
Li(n£) «-> Li + (ls 2 ) + H is treated accord- 
ing to Barklem, et al. (2003). More details 
are given in Sbordone et al. (2010). Finally, 
3D non-LTE synthetic line profiles of the Li i 
A 670.8 nm doublet are computed with the line 
formation code Linfor3D 2 , using the depar- 
ture coefficients b t = n,(NLTE)/n,(LTE) pro- 
vided by NLTE3D for each level / of the 
lithium model atom as a function of geometri- 
cal position within the 3D model atmospheres. 
At this stage, the 3D hydrodynamical velocity 
field provides the differential line shifts along 
each line of sight. As demonstrated in Steffen 
et al. (2010b), non-LTE effects are very im- 
portant for the 3D models of the metal-poor 
dwarfs considered here: they strongly reduce 
the height range of line formation such that 
the 3D non-LTE equivalent width is smaller 
by roughly a factor 2 compared to 3D LTE. 
Ironically, the line strength predicted by stan- 
dard ID mixing-length models in LTE are 
close to the results obtained from elaborate 3D 
non-LTE calculations. Nevertheless, when it 



2 http: //www. aip.de/~mst/linfor3D_ 
main.html 



comes to the determination of the 6 Li/ 7 Li iso- 
topic ratio, it is important to account for the 
convective line asymmetry, which requires the 
full 3D non-LTE line formation calculations. 

3. Theoretical 3D non-LTE correction 
of 1D LTE 6 Li/ 7 Li determinations 

As outlined above, the 6 Li abundance is sys- 
tematically overestimated if the intrinsic asym- 
metry of the 7 Li line components is ignored. 
To quantify this bias theoretically, we rely on 
synthetic spectra. The idea is as follows: we 
represent the observation by the synthetic 3D 
non-LTE line profile of the 7 Li line blend, com- 
puted with zero 6 Li content. The non-thermal 
line broadening is naturally provided by the 3D 
hydrodynamical velocity field, which replaces 
the classical concept of micro- and macro- 
turbulence, and also gives rise to a convec- 
tive blue-shift and an intrinsic line asymme- 
try. Optionally, the line blend is subsequently 
broadened by rotation (u sin z) and a Gaussian 
instrumental profile (FWHM denoting its full 
width at half maximum). Now this realistic 
proxy of the 7 Li line blend is fitted by a grid of 
'classical' ID LTE synthetic line profiles. The 
ID LTE line profiles are computed from so- 
called LHD models, ID mixing-length model 
atmospheres that have the same stellar param- 
eters and employ the same microphysics and 
radiative transfer scheme as the corresponding 
3D model, assuming a mixing-length param- 
eter of aMLT = 0.5, and a depth-independent 
microturbulence of ^ m i c = 1.5km/s. 

Five parameters are varied independently 
until the best fit is found (x 2 minimization): 
in addition to the total 6 Li+ 7 Li abundance, 
A(Li), and the isotopic ratio, <7(Li)= 6 Li/ 7 Li, 
which control line strength and line asymme- 
try, respectively, we also allow for a residual 
line broadening described by a Gaussian ker- 
nel with half-width Vbr, a global line shift, Ay, 
and a global scaling of the continuum level, S c . 
For technical reasons, the rotational line broad- 
ening is not treated as a fitting parameter; it is 
fixed to the value used in the 3D spectrum syn- 
thesis. 

Note that, although the five fitting param- 
eters are correlated, each one has a distinctly 
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Fig. 1. Example of a 3D non-LTE synthetic profile of the Lii 670.8 nm feature, for stellar parameters 
r eft = 6300 K, \ogg = 4.0, [Fe/H]=-2. Increasing the 6 Li content from to 5% leads to an extra depression 
of the red wing by about 0.5% (indicated by the gray area), leaving the blue wing unaffected. In order to 
detect such a small effect, observed spectra are required to have an S/N in excess of ss 400. 



different effect on the line profile. The broaden- 
ing Vbr is fixed by the slope of the blue wing, 
which is not affected by the 6 Li content, (see 
Fig. 1), while q(Li) is fixed by the shape of the 
red wing. This allows to obtain unambiguously 
a unique solution although a given half width 
of the line blend can be obtained by different 
combinations of q(Li) and Vbr, but with higher 
X 2 values than the best solution. 

Finally, the 3D non-LTE correction of the 
6 Li/ 7 Li isotopic ratio is defined as the differ- 
ence between the values of q{lj) that provide 
the best fit (*) to the given profile using ID 
LTE and 3D non-LTE profiles, respectively: 
A^(Li) = 4* 1DLTE) (Li) - 4 ( * 3DNLTE) (Li) 3 . This 
differential correction is meant to be subtracted 
from the 6 Li/ 7 Li isotopic ratio derived from the 
ID LTE analysis to correct for the bias intro- 
duced by neglecting the intrinsic line asym- 
metry: <?(3DNLTE)(Li) = <?(lDLTE)(Li) - Atffji). 



9(3 DNLT m(Li) is zero by construction here 



Note that this correction procedure properly 
accounts for radiative transfer in the lines, in- 
cluding saturation effects. 

We have determined Ag(Li) according to 
the method outlined above for each of the 
3D model atmospheres in our grid, consid- 
ering two cases: A^(A>(Li) and A^ ( B)(Li) are 
the result of fitting the unprocessed (v sin i = 
0.0 km/s, FWHM = 0.0 km/s) and the broad- 
ened (v sin / = 1.5 km/s, FWHM = 2.5 km/s) 
synthetic line profile, respectively. Case (B) is 
the more realistic one, while case (A) has been 
included only to show how the amount of in- 
strumental broadening is affecting the ampli- 
tude of Ag(Li). Note also that the 3D correc- 
tions given in previous studies by (Steffen et 
al. 2010a,b) correspond to case (A). 

The results of this computationally very ex- 
pensive procedure are given in cols. (8) and 
(9) of Table 1 . They may be aproximated fairly 
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Table 1. List of 3D models used in the present study. Columns (2)-(6) give effective temper- 
ature, surface gravity, metallicity, number of opacity bins used in the radiation hydrodynamics 
simulation, and number of snapshots selected for spectrum synthesis. The equivalent width of 
the synthetic 3D non-LTE 7 Li doublet at A 670.8 nm, assuming A(Li)=2.2 and no 6 Li, is given 
in col. (7). Columns (8) and (9) show Ag(A)(Li) and Aq^ilA), the 3D non-LTE correction of the 
6 Li/ 7 Li isotopic ratio inferred from fitting the unprocessed (no rotational and instrumental broad- 
ening) and the broadened (v sin/ = 1.5 km/s, FWHM = 2.5 km/s) 3D non-LTE line profile, 
respectively, with a grid of ID LTE profiles. 



Model 


f l) 
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[K] 


log g 
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Nbin 


Nsnap 


W 2 ' 
[mA] 


A^(Li) 
(A) 


^(IDLTE)^ 1 ) 
[X100] 

(B) 


d3t59g35mm30n01 


5873 


3.5 


-3.0 


6 


20 


41.7 


2.13 


1.99 


d3t59g40mm30n02 


5846 


4.0 


-3.0 


6 


20 


44.9 


0.88 


0.86 


d3t59g45mm30n01 


5924 


4.5 


-3.0 


6 


19 


39.9 


0.63 


0.61 


d3t63g35mm30n01 


6306 


3.5 


-3.0 


6 


20 


21.7 


3.44 


3.12 


d3t63g40mm30n01 


6269 


4.0 


-3.0 


6 


20 


23.9 


1.51 


1.38 


d3t63g45mm30n01 
d3t65g40mm30n01 


6272 
6408 


4.5 
4.0 


-3.0 
-3.0 


6 
6 


18 

20 


24.3 
20.0 


0.85 
1.70 


0.80 
1.56 


d3t65g45mm30n01 


6556 


4.5 


-3.0 


6 


12 


16.4 


1.25 


1.16 


d3t59g35mm20n01 


5861 


3.5 


-2.0 


6 


20 


42.8 


1.88 


1.76 


d3t59g40mm20n02 


5856 


4.0 


-2.0 


6 


20 


45.2 


0.73 


0.69 


d3t59g45mm20n01 


5923 


4.5 


-2.0 


6 


18 


42.3 


0.45 


0.41 


d3t63g35mm20n01 


6287 


3.5 


-2.0 


6 


20 


22.1 


4.04 


3.75 


d3t63g40mm20n01 


6278 


4.0 


-2.0 


6 


16 


23.7 


1.48 


1.35 


d3t63g45mm20n01 


6323 


4.5 


-2.0 


6 


19 


23.0 


0.94 


0.88 


d3t65g40mm20n01 
d3t65g45mm20n01 


6534 
6533 


4.0 
4.5 


-2.0 
-2.0 


6 
6 


19 
19 


16.2 
17.0 


2.22 
1.21 


2.01 
1.12 


d3t59g35mml0n02 


5890 


3.5 


-1.0 


6 


20 


38.0 


2.91 


2.71 


d3t59g40mml0n02 


5850 


4.0 


-1.0 


6 


20 


41.5 


1.45 


1.36 


d3t59g45mml0n01 


5923 


4.5 


-1.0 


6 


08 


38.2 


0.83 


0.76 


d3t63g35mml0n01 
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3.5 


-1.0 


6 


20 


23.0 


4.57 


4.18 


d3t63g40mml0n01 


6261 


4.0 


-1.0 


6 


20 


22.0 


2.33 


2.15 


d3t63g45mml0n01 


6238 


4.5 


-1.0 


6 


20 


23.4 


1.23 


1.14 


d3t65g40mml0n01 


6503 


4.0 


-1.0 


6 


20 


15.5 


3.14 


2.82 


d3t65g45mml0n01 


6456 


4.5 


-1.0 


6 


19 


17.1 


1.44 


1.33 


d3t59g35mm00n01 
d3t59g40mm00n01 


5884 
5928 


3.5 
4.0 


0.0 
0.0 


5 
5 


20 
18 


33.0 
31.1 


3.37 
1.83 


3.06 
1.70 


d3t59g45mm00n01 


5865 


4.5 


0.0 


5 


19 


34.7 


1.15 


1.06 


d3t63g40mm00n01 


6229 


4.0 


0.0 


5 


20 


19.1 


2.66 


2.43 


d3t63g45mm00n01 


6233 


4.5 


0.0 


5 


20 


19.3 


1.49 


1.47 


d3t65g40mm00n01 


6484 


4.0 


0.0 


5 


20 


12.8 


3.88 


3.55 


d3t65g45mm00n01 


6456 


4.5 


0.0 


5 


20 


13.8 


1.75 


1.61 



Notes: " averaged over selected snapshots; 2) averaged over selected 3D non-LTE spectra 



well by the following polyonmial expression: 

A4( B )(Li){7/ eff , log g, [Fe/H]} = a m + a 01 Z + 
(flio + a u Z) Y + (fl 2 o + «2i Z) Y 2 + 
(a 30 + a 3i Z) X + (a 40 + a 4 i Z) X Y + 



(a 5 o + a 5l Z + a 52 Z 2 ) X 2 (1) 

where X=(T eS - 6300)/1000, Yslogg - 4.0, 
and Zs[Fe/H]+2. The coefficients a (J - are listed 
in Table 2. Contour plots of this polynomial ap- 
proximation are shown in Fig. 2. 
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Fig. 2. Contours of Ag (B )(Li) in the T e g - log g plane at constant metallicity [Fe/H] = -3, -2, -1, and 0, 
respectively (from top left to bottom right) according to the approximation given by Eq. (1) and Table 2. 
White symbols mark the positions of the stars from the A06 sample with Fe/H] < -2.5 (top left), -2.5 < 
[Fe/H] < -1.5 (top right), and -1.5 < [Fe/H] < -0.5 (bottom left). 

Table 2. Coefficients representing A^(B)(Li) The results for [Fe/H]=0.0 are most uncertain, 
as a function of T e s, log g, and [Fe/H] accord- since they are based on fewer data points, 
ing to Eq. (1). 



j 


a 0j x 100 


ay x 100 


a 2j x 100 





1.6984 


0.4311 


0.0000 


1 


-2.6658 


-0.4547 


0.0000 


2 


2.3370 


-0.3211 


0.0000 


3 


2.3490 


0.6300 


0.0000 


4 


-2.7366 


-0.6089 


0.0000 


5 


0.1786 


0.5657 


0.5385 



At given metallicity, the corrections are 
largest for low gravity and high effective tem- 
perature. They increase towards higher metal- 
licity. Note that the corrections are strictly 
valid only for a Li abundance of A(Li)=2.2. 





0.00 
6Li/7Li 



Fig. 3. Histogram of the 6 Li/ 7 Li determinations by 
A06 before (top) and after (bottom) application of 
the 3D correction Ag ( B)(Li). 



The analysis of A06 utilizes ID LTE pro- 
files computed from MARCS model atmo- 
spheres. Hence, the correction Ag(£)(Li) should 
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be applied to their 6 Li/ 7 Li isotopic ratios. The 
resulting downward corrections are typically in 
the range 1% < Ag(£)(Li) < 2% for the stars 
of their sample (cf. Fig. 2). After subtracting 
the individual Ag (fi) (Li) for each of these stars, 
interpolated via Eq. (1) to their T e g, log g, and 
[Fe/H], the mean 6 Li/ 7 Li isotopic ratio of the 
sample is reduced from 0.0212 to 0.0088. The 
center of the distribution of 6 Li/ 7 Li is shifted 
essentially to zero (see Fig. 3). 

Based on the original 6 Li/ 7 Li isotopic ra- 
tios and their l<x error bars as determined by 
A06, the number of stars with a 6 Li detection 
above the 2<x and 3cr level is 9 and 5, respec- 
tively, out of 24. After correction, the number 
of 2cr and 3<x detections is reduced to 5 and 2, 
respectively. The 3<x detections after correction 
are HD 106038 and G020-024, the remaining 
2cr detections are HD 102200, HD 160617, and 
CD-30 18140. 

We note that HD 106038 survives as a 3<x 
detection only because of its amazingly small 
error bar of +0.006. Assuming a more typi- 
cal error of +0.013, this object would not even 
qualify as a 2<x detection. HD 102200 seems 
to be a clear 2<x detection (<7(Li) = 0.036 + 
0.013), while CD-30 18140 is a marginal case 
(q(U) = 0.029 ± 0.013). The remaining ob- 
jects, G020-024 and HD 160617, are discussed 
in more detail in Sect. 5, where their spectra 
are reanalyzed with 3D non-LTE line profiles; 
none of them provides convincing evidence for 
the presence of 6 Li. 

4. Fixing the residual line broadening 
by 'calibration lines' 

So far we have derived all information from fit- 
ting the Li 670.8 nm feature, for which we have 
a 3D non-LTE line formation model. In princi- 
ple, the accuracy of the fitting procedure can 
be improved by fixing the value of the resid- 
ual line broadening, Vbr, from additional un- 
blended 'calibration lines'. This approach has 
been adopted by A06 (see also Lind et al. 
2012, henceforth L12). This fact complicates 
the comparison with our results, and the cor- 
rection procedure described in Sect. 3 explains 
only part of the differences. We show in the fol- 
lowing that 'calibration lines' may introduce 



Table 3. List of calibration lines. Columns (2) 
to (4) show their wavelength, excitation po- 
tential, and equivalent width. The residual line 
broadening derived from ID LTE fitting and 
the resulting error of 6 Li/ 7 Li is given in cols. 
(5) and (6). The last row refers to the Li line; 
cols. (5) and (6) show the results of 3D non- 
LTE fitting without calibration lines. 



Ion 


A 








6q (W > 




[nm] 


[eV] 


[pm] 


[kms] 


[X100] 


Fei 


538.337 


4.312 


2.74 


4.88 


0.35 


Fei 


606.548 


2.608 


1.38 


4.63 


1.25 


Fei 


613.662 


2.453 


2.47 


4.52 


1.64 


Fei 


613.769 


2.588 


1.87 


4.59 


1.39 


Fei 


623.072 


2.559 


2.51 


4.54 


1.57 


Fei 


625.256 


2.404 


1.64 


4.56 


1.50 


Fei 


639.360 


2.433 


2.53 


4.51 


1.67 


Fei 


649.498 


2.404 


3.54 


4.39 


2.08 


Fei 


667.799 


2.692 


1.58 


4.62 


1.29 


Cai 


671.768 


2.709 


0.78 


4.83 


0.53 


Li i 


670.784 


0.000 


2.37 


4.97 


0.00 



additional uncertainties, unless their 3D non- 
LTE line formation is fully understood. 

In the framework of ID models, V B r, rep- 
resents the combined effect of macroturbulence 
and instrumental broadening (for given micro- 
turbulence and rotational velocity) and is ex- 
pected to be independent of the spectral line. In 
general, an average value of Vbr from several 
calibration lines is then used for the analysis of 
the Li line, thus reducing the number of free 
fitting parameters, and hence the formal errors 
of the fitting results. 

We have simulated this procedure using 
synthetic lines only. Our selection of 10 cali- 
bration lines is compiled in Tab. 3. The lines 
have been selected to have similar wave- 
length, equivalent width, and excitation en- 
ergy (relative to the ionization continuum) as 
Li 1 670.8 nm. Synthetic line profiles generated 
from 3D model d3t63g40mm20n01 (assuming 
LTE) represent the observed line profiles and 
are fitted with ID LTE profiles from the cor- 
responding LHD model to derive V B r. The re- 
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Table 4. List of observed Li I A 670.8 nm spectra of eight metal-poor turn-off stars, analyzed in 
this work for the presence of 6 Li. 



Star 


T eB [K] 


log 9 


[Fe/H] 


R=A/AA 


S/N "> 


Instrument 


Reference 


HD 74000 


6203 


4.03 


-2.05 


120000 


600 


ES03.6 / HARPS 


C07 *> 


G271-162 


6330 


4.00 


-2.25 


110000 


600 


VLT / UVES 


N00 2) 


HD 84937 


6310 


4.10 


-2.40 


100000 


650 


CFHT/ GECKO 


C99 3) 


G020-024 


6247 


3.98 


-1.89 


120000 


420 


VLT / UVES 


A06 4) 


HD 140283 


5753 


3.70 


-2.40 


95 000 


1000 


SUBARU /HDS 


A04 5) 


HD 160617 


5990 


3.79 


-1.76 


100000 


400 


ES03.6 / HARPS 


M12 6) 


G64-12 


6371 


4.26 


-3.24 


95 000 


620 


Keck/ HIRES 


A08 7) 


G275-4 


6338 


4.32 


-3.21 


93 000 


700 


VLT / UVES 


N09 8) 



Notes: 11 Cayrel et al. (2007), 21 Nissen et al. (2000), 3) Cayrel et al. (1999), 4) Asplund et al. (2006), 51 Aoki et al. (2004), 
f,) L. Monaco and G. Lo Curto (2012, priv. comm.), 7) Asplund & Melendez (2008), 8) P.E. Nissen (2009, priv. comm.), 
*' S/N estimated from the (re-reduced) spectra actually used in this work, may differ from literature values. 



suits are given in col. (5) of Tab. 3. These num- 
bers have to be compared with the value of 
Vbr obtained by fitting the 3D non-LTE Li pro- 
file of the same 3D model with ID LTE pro- 
files from LHD models, V BR = 4.97 km/s (last 
row of Tab. 3). When fixing Vbr to the val- 
ues implied by the individual calibration lines, 
the resulting 6 Li/ 7 Li isotopic ratio changes by 
&7 (1D) (Li) with respect to the value obtained 
when treating V B r as a free fitting parameter. 

We find that, even under these idealized 
conditions, 6q (lD \Li) is not constant but varies 
from line to line (see col. (6) of Tab. 3). This is 
an indication that ID micro / macro model is 
not perfectly adequate to describe the 3D hy- 
drodynamical velocity field. Moreover, since 
Sq^ D \lA) is positive for all of our calibration 
lines, the use of calibration lines leads to an ad- 
ditional overestimation of the Li isotopic ratio 
by up to 2%. The total error of the ID analy- 
sis is obtained by adding the corrections from 
Tables 1 and 3, Ag(Li) = Aq (B )(U) + Sq (W) (U). 

In a further experiment, we have fitted both 
the calibration lines and the Li feature with 3D 
LTE line profiles. By construction, the derived 
Vbr is then equal to FWHM of the assumed 
instrumental profile, and is identical for all cal- 
ibration lines. The Li isotopic ratio obtained 
from fitting the 3D non-LTE Li feature with 
3D LTE profiles and fixed Vbr = FWHM is 
&7 (3D) (Li)=0.054, compared to 0.003 if V BR is 
a free fitting parameter. The basic reason for 
this alarming result is that the half width of 



the Li line is significantly smaller in LTE than 
in non-LTE, which must be compensated by 
a higher 6 Li content. This effect may explain 
why A06 find higher 6 Li/ 7 Li isotopic ratios in 
their 3D LTE analysis (with respect to ID LTE, 
their Table 5). 

We note that the above results are con- 
sistent with previous findings by Steffen et 
al. (2010b), who analyzed the spectrum of 
HD 74000 with a subset of 6 clean Fei cali- 
bration lines, and obtained significantly higher 
6 Li/ 7 Li isotopic ratios compared to the analysis 
without calibration lines. The same behavior 
is seen by LI 2, if they fit both the calibration 
lines and Li with 3D LTE profiles. Using in- 
stead non-LTE profiles for fitting both Ca i and 
Li i lines leads to consistent 6 Li/ 7 Li isotopic ra- 
tios for HD 84937. 



5. Analysis of observed spectra: 3D 
non-LTE versus 1D LTE line fitting 

With the results given in Sect. 3, the recom- 
mended method of analysis of the Li 1 670.8 nm 
doublet in an observed spectrum is as follows. 
The first step is to produce a grid of ID LTE 
synthetic line profiles as a function of A(Li), 
q(Li), using any standard ID mixing-length 
model with the correct stellar parameters (r e ff, 
log g, [Fe/H]). This grid is then used to find the 
best fit to the observed line profile, as described 
above (no calibration lines). Finally, the result- 
ing q(Li) is corrected for 3D effects by sub- 
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HD 84937 



1.02 



1.00 



£ 0.98 



E 

o 0.96 



0.94 



0.92 



Fitting by: d3t63g40mm20n01-16ss-3D (non-LTE) 



A(Li) 
" 6 Li/ 7 Li 
shift 
xi 



2.20 

0.0512 +/-0.022& 
-0.0217 km/s \ 
1 .86 +/f 0.40 km/s 



" fwhM 


= 0.00 krti/s 


v sinii | 


= 1 .00 kfn/s 


_ clevel \ 


= 0.999537 


_ chi^sq ' 


^ 34.8 ' 


-niter 




- npixel 


' :V'¥ 
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" int2 


2 



-40 
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0.98 
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0.90 
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Fig. 4. The best fitting 3D non-LTE Li 670.8 nm line profile (green) superimposed on the observed spec- 
trum (solid black) of HD 84937 (top) and G020-024 (bottom). The difference between computed and ob- 
served spectrum (xlO; dashed line) indicates the quality of the fit, which is much better for HD 84937 than 
for G020-024, the latter suffering from a 'perturbation' in its lower red wing. 



tracting the differential correction Aq^Li) de- 
fined in Sect. 3, interpolated to the actual stellar 
parameters according to Eq. (1) and Table 2. 

In the following, however, we disregard 
this differential approach, and present instead 
preliminary results of fitting the observed Li i 
A 670.8 nm spectra of eight halo turn-off stars 
(see Table 4) directly with our grids of 3D non- 
LTE and ID LTE synthetic line profiles, 
computed from 3D CO 5 BOLD models and 
ID LHD models, respectively. As described 



above, five free fitting parameters (A(Li), ^(Li), 
Vbr, Ay, S c ) are varied independently to find 
the best fit. Formally, we fix the rotational ve- 
locity to v sin i— 1.0 km/s, noting that the de- 
rived 6 Li/ 7 Li isotopic ratio is insensitive to this 
assumption (see Steffen et al. 2010b). 

The results are presented in Table 5. The er- 
rors quoted for Vbr and g(Li) are the formal 
1 <x confidence intervals due to the finite S/N 
of the fitted spectra. They are computed in the 
usual way (see e.g. Nissen et al. 2000; Asplund 
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Table 5. Fitting the observed spectra of Table 4 with 3D non-LTE and ID LTE line profiles. 



Star Model synthetic A (Li) » q(Li) [xlOO] 

name spectrum [km/s] this work A06,A08,L12 



HD 74000 d3t63 g 40mm20n01 3D NLTE 2.25 2.9 ±0.4 -0.8 ±1.2 
t6280 g 40mm20a05ob6 ID LTE 2.23 5.7 ± 0.2 0.6 ±1.1 



(jZ / 1 — 1 oz 


d3t63g40mm20n01 


1T\ Ml TTT 


L.Ly 


j.J ± U.J 


U.o ±1.1 






t6280g40mm20a05ob6 


1 T^ T TT-T 
1JJ Llll 


2. LI 


a 1 _i_ n i 
O.l ± U.2 


O 1 _i_ 1 

2.1 ± 1.2 


i.y ± 1.2 ; 


HD 84937 


d3t63g40mm20n01 


3D NLTE 


2.20 


3.1 ±0.7 


5.1 ±2.3 


0.7 ± 0.5 5) 




t6280g40mm20a05ob6 


ID Lit, 


L. to 


J.o ± U.4 


6.6 ± LA 




G020-024 


d3t63g40mm20n01 


3D NLTE 


2.21 


0.9 ± 1.3 


9.9 ± 2.4 


— 




t6280g40mm20a05ob6 


ID LTE 


2.19 


5.0 ±0.4 


11.7 ±2.5 


7.0 ± 1.7 3) 


HD 140283 


d3t57g37mm20n01 


3D NLTE 


2.20 


3.7 ±0.1 


-0.6 ± 0.4 


-0.4 ± 0.2 5) 




d3t57g37mm30n01 


3D NLTE 


2.22 


3.3 ±0.1 


-1.0 ±0.4 






t5775g37mm20ml3a05ob6 


ID LTE 


2.18 


5.5 ±0.1 


0.5 ± 0.4 


0.8 ± 0.6 3) 




t5789 g 37mm30ml3a05ob6 


ID LTE 


2.19 


5.5 ±0.1 


0.4 ± 0.4 




HD 160617 


d3t59g40mm20n02 


3D NLTE 


2.14 


4.5 ± 0.2 


-0.5 ± 1.2 






d3t59g35mm20n01 


3D NLTE 


2.17 


3.3 ±0.3 


-1.3 ± 1.2 






t5860g40mm20ml3a05 


ID LTE 


2.11 


5.8 ±0.2 


0.1 ± 1.2 


3.6 ± 1.0 3) 




t5860g35mm20a05 


ID LTE 


2.13 


5.8 ±0.2 


0.3 ± 1.2 




G64-12 


d3t63g40mm30n01 


3D NLTE 


2.19 


5.0 ±0.5 


0.8 ± 2.7 


2.1 ± 1.2 5) 




d3t63g45mm30n01 


3D NLTE 


2.19 


5.5 ± 0.5 


1.2 ±2.8 






t6270 g 40mm30ml3a05 


ID LTE 


2.16 


7.1 ±0.4 


2.3 ± 2.9 


5.9±2.1 4) 




t6270 g 45mm30ml3a05ob6 


ID LTE 


2.14 


7.0 ± 0.4 


2.0 ± 2.9 




G275-4 


d3t63g40mm30n01 


3D NLTE 


2.11 


3.6 ±0.6 


3.5 ± 2.4 






d3t63g45mm30n01 


3D NLTE 


2.11 


4.3 ± 0.5 


3.8 ± 2.4 






t6270g40mm30ml3a05 


ID LTE 


2.08 


6.3 ±0.4 


4.7 ± 2.5 






t6270 g 45mm30ml3a05ob6 


ID LTE 


2.06 


6.2 ± 0.4 


4.4 ± 2.4 





Notes: " log [n( f 'Li) + n( 7 Li)] - log «(H) + 12; 2) Gaussian kernel, 3 » Asplund et al. (2006), 4) Asplund & Melendez (2008), 
5) Lind etal. (2012) 



et al. 2006) by finding the distance Ap of the 
parameter of interest p from is optimum value 
p Q such that^ 2 (/? + Ap) -x 2 (po) = Ay 2 = 1, 
fixing p and minimizing Ay 2 over the remain- 
ing free fittting parameters. Here x 2 is defined 
as 

i 

where (9, and 5, are the observed and synthetic 
flux at wavelength point i, and cr=(S/N)^ 1 , with 
S/N taken from Tab. 4. Note that this is a lower 
limit to the real error, which may have several 
other significant contributions (see discussion 
in Garcia Perez et al. 2009). 



As expected, the 3D analysis yields sys- 
tematically lower 6 Li/ 7 Li isotopic ratios by up 
to 1.8%. For comparison, we list in the last 
column the ID LTE results of A06, A08, and 
the 3D non-LTE results of L12 from their case 
(a) for the stars in common with our sam- 
ple. In some cases, the agreement is very good 
(G27 1-162, HD 140283), even though the anal- 
ysis is based on different observational data. 
In other cases, we obtain significantly lower 
6 Li/ 7 Li ratios (HD 160617, G64-12). The most 
striking discrepancy, this time in the opposite 
direction, is seen for HD 84937. 

G020-024 is a special case. We have re- 
trieved the spectra obtained by A06 from the 
ESO archive and re-reduced them with the cur- 
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rent UVES pipeline. We note that all six sub- 
exposures show an unexplained 'perturbation' 
in the lower red wing of the Li 670.8 nm line 
that cannot be fitted properly, neither with ID 
LTE nor with 3D non-LTE model spectra. The 
best 3D non-LTE fit suggests a very high 6 Li 
content near 10% (see Fig. 4), even exceeding 
the ID LTE result by A06 of 7%. In view of 
the poor quality of the fit, we consider this re- 
sult as doubtful. It is worth mentioning that 
G020-024 is listed as a suspected binary in 
Fouts (1987); moreover, A06 report that this 
star shows an unusually large discrepancy be- 
tween the photometric and the H a -based T e g. 
Possibly, G020-024 is a spectroscopic binary 
whose components are of similar spectral type. 

Disregarding G020-024, we are left with 
one formal 2cr detection (HD 84937), 
one 1 <x detection (G275^), and five 
non-detections (HD 74000, G271-162, 
HD 140283, HD160617, G64-12) out of the 
eight stars, when considering the 3D non-LTE 
results only. In ID LTE, G271-162 would turn 
into a formal 1 cr detection (see Tab. 5). 

6. Conclusions 

The 6 Li/ 7 Li isotopic ratio derived by fitting the 
Li I doublet with 3D non-LTE synthetic line 
profiles is shown to be 1% to 2% lower than 
what is obtained with ID LTE profiles. This 
result implies that only 2 out of the 24 stars 
of the Asplund et al. (2006) sample would for- 
mally remain significant (3 cr) 6 Li detections 
when subjected to a 3D non-LTE analysis. 

In another theoretical case study we have 
demonstrated that the difference between 3D 
non-LTE and ID LTE results increases even 
more if we would rely on additional 'cali- 
bration lines' to fix the residual line broad- 
ening, as advocated e.g. by A06, A08, L12. 
The number of possible 6 Li detections is thus 
further reduced. We conclude that the usage 
of additional 'calibration lines', even if care- 
fully selected, introduces additional uncertain- 
ties rather than reducing the error of the anal- 
ysis, unless the 3D non-LTE line formation is 
fully understood for all involved lines. 

Finally, we have analyzed available high 
quality spectra of eight turn-off halo stars, both 



with ID LTE and 3D non-LTE modeling. In 
most cases, the evidence for the presence of 6 Li 
is not significant. Only in the case of HD 84937 
it is difficult to deny the signature of Li. 
Surprisingly, L12 no longer find any signifi- 
cant evidence for the presence of 6 Li in this 
object, which since almost two decades has 
been considered as an undisputed 6 Li detec- 
tion. From the results reported by LI 2, it is 
unclear whether the reason for this disturbing 
discrepancy is related to the new superior ob- 
servational data or to the improved method of 
analysis. 

We conclude that 3D model atmospheres 
can indeed help to solve the second lithium 
problem. In view of the 3D non-LTE results 
reported in this work (and by LI 2), it seems 
that the presence of 6 Li in the atmospheres of 
galactic halo stars is rather the exception than 
the rule, and hence does not necessarily consti- 
tute a cosmological problem. 
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